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Abstract 

The integrated luminosity ~ 4000n6"^ of around </> meson mass (5.0 x 10^ 
of </>'s) has been cohected with the CMD-2 detector at the VEPP-2M cohider. 
A latest analysis of the KsKl coupled decays based on 30% of available data is 
presented in this paper. 

The KsKl pairs from (p decays were reconstructed in the drift chamber when 
both kaons decayed into two charged particles. From a sample of 1423 coupled 
decays a selection of candidates to the CP violating K^ -^ tt+tt" decay was 
performed. CP violating decays were not identified because of the domination of 
events with a K^ regenerating at the Be beam pipe into Ks and a background 
from Kl semileptonic decays. 

The regeneration cross section of 110 MeV/c K^s was found to be a^fg = 
53 lb 17 mb in agreement with theoretical expectations. The angular distribution 
of Ks^s after regeneration and the total cross section of K^ for Be have been 
measured. 



1 Introduction 



As was realized at the very early steps of the (j) meson studies at the coUiding beam 
machines, KsKl pairs ( ~ 34% of all decays) can be used as a source for studying 
CP and CPT violation. These suggestions, including studies of quantum mechanical 
correlations, were discussed for experiments at the Novosibirsk electron-positron col- 
lider VEPP-2M |I|-[^. The coupled decays of the KsK^ pairs will allow demonstration 
of the quantum mechanical correlations of the two particle decays (Einstein-Podolsky- 
Rosen paradox) @] . 

The (j) resonance produced in e"*"e~ collisions is also a source of low momentum 
neutral and charged kaon pairs not available from other sources for studies of nuclear 
interactions. Mesons in each pair are produced with opposite and equal momenta 
and a detector with good resolution and reconstruction efficiency allows one to use one 
reconstructed kaon as a tag for another. The most natural way is to use a reconstructed 
decay of Ks — > tt+tt^ as a tag for K^. In this case the momentum and direction of K^ 
are completely determined. 

The idea of creating an intensive source of (f) mesons has been discussed by many 
authors P, |§|. The flux of events at these so-called " 0-factories" now under construction 
1^, Q will make feasible new precise measurements of a possible direct component in 
the decay Kl -^ tt+tt", tt^tt" {^'/^), as well as an observation of CP- violating three 
pion decays of the Ks- Studies of the oscillations in the joint decay distributions could 
provide information about real and imaginary parts of any CP-violating amplitude. 

At the VEPP-2M collider at Novosibirsk, which could be considered as a pre 0- 
factory, we have been running with the CMD-2 detector preparing for experiments at 



the 0-factory which is under construction here. Studies of an upgraded detector and ac- 
celerator are in progress, including an intermediate 10^^cr7i~^s~^ luminosity collider for 
investigating the idea of the round beams, an important ingredient of the Novosibirsk 
0- factory project [^ |T0|. 

2 The CMD-2 Detector 



The CMD-2 detector has been described in more detail elsewhere y, [Tl[ . The main 
systems of the detector are shown in Figure 1. 

The Csl barrel calorimeter with a 6 x 6 x 15 cm'^ crystal size is placed outside of a 
0.4 r.l. superconducting solenoid with a 1 Tesla magnetic field. The endcap calorimeter 
is made of 2.5 x 2.5 x 15 cm^ BGO crystals and has not been installed for the data 
presented here. 

The drift chamber(DC) with a 30 cm outer radius and a 44 cm length has about 
250 /im resolution transverse to the beam and 0.5 cm longitudinally and is placed 
inside the solenoid. The vertex reconstruction resolution for the neutral kaon decays 
into charged particles is about 0.15 cm radially. The muon range system consists of 
streamer tubes and has 1-3 cm spatial resolution. 

A 3.4 cm diameter vacuum beam pipe is made of Be with a 0.077 cm wall thickness 
and may be considered as a target for studies of the kaon nuclear interaction. 

A data sample of the integrated luminosity of 1500 nb~^ has been analyzed around 
the corresponding to about 1.7 x 10^ produced 0's. About 300 nb~^ were used to 
measure the meson parameters and its branching ratios into four major decay modes 
[ 1^ . The largest part of the integrated luminosity ~ 1200 nb~^ was used for the studies 



of rare decay modes of 0. Some preliminary results were published in 

These data were also intensively used for studies of the detector performance and 
software development. Several versions of the reconstruction program have been tested 
and a reconstruction efficiency of about 95-97% per charged track has been achieved 
with a momentum uncertainty of 4% for 250 MeV/c charged particles. The energy 
resolution of 8% for photons in the Csl calorimeter has been obtained. 

The number of reconstructed KsKl coupled decay events is now twice the sample 
in the preliminary publication [jl3|. The latest analysis of the KsK^ coupled decays is 
presented in this paper. 

3 Selection of KsKi coupled decays 

Candidates were selected from a sample in which two vertices, each with two 
oppositely charged tracks, were observed within 15 cm from the beam axes and all 





Figure 1: Horizontal and vertical cross sections of the CMD-2 detector. 1 - vacuum 
chamber; 2 - drift chamber; 3 - Z-chamber; 4 - main solenoid; 5 - compensating solenoid; 
6 - BGO endcap calorimeter; 7 - Csl barrel calorimeter; 8 - muon range system; 9 - 
magnet yoke; 10 - collider lenses. 



tracks were reconstructed. An example of such an event is shown in Fig. ^. Figure |^a 
shows a scatterplot of the invariant mass Minv of two charged tracks assuming they are 
pions vs. the missing momentum Pmis (which is equivalent to kaon momentum P^ in 
case of two pion decay) for the particle which decays first. Figure ^d shows the same 
plot for a particle from the second vertex. The concentration corresponding to Ks 
mesons decaying into a pair of charged pions dominates in the Fig. ^ and is seen in 
Fig. Pd showing events in which Ks appears in the second vertex. The cuts 470 MeV 
< Minv < 525 MeV and 80 MeV/c < Pr < 140 MeV/c with an additional requirement 
to have another reconstructed vertex in the Pmis direction select events with Ks in one 
of the vertices. The K^ is expected to be in the other one. The above cuts represent 
three standard deviations of the detector resolution. 

Figures |^ and ^jd show peaked Minv and Pk distributions corresponding to Ks 
after above cuts (shaded). The Mj„j, and Pk parameters for the other vertex give broad 
distributions expected for the main 3-body Kl decays (into vr/iz/, vrez/, vt+tt^tt^). 

Figure 4a shows the decay length distribution for selected Ks- The exponential 
decay length is seen with a correct value 0.55±0.02 cm and with a vertex position 
resolution of 0.15±0.03 cm. 

The total number of the reconstructed K^ vertices accompanying Ks decays was 
found to be 1423. The interpretation of the decay radius distribution for the Kl shown 





Figure 2: ^ KsKl event with a coupled decay. 



in Figure 4b is more complicated. 

The decay length for Kl is about 350 cm, so only about 4% of KsKls both decay 
within 15 cm of DC and only half of them both decay into charged modes. With the DC 
angular acceptance one can expect to detect 1% of KsKi pairs if the reconstruction 
efficiency is radially uniform. 

But the reconstruction efficiency in the jet type DC cells drops down for tracks with 
a vertex radius greater than approximately 5 cm leaving only 0.64% of all detected 
events with Ks and Kl vertices. To describe the Kl decay radius distribution, it is 
assumed that the reconstruction efficiency is uniform to a certain radius and then drops 
down exponentially. 

A clear peak at a radius of 1.7 cm corresponds to the Kls which interacted with 
nuclei in the Be tube. This peak has two contributions - Kl regeneration into Ks with 
Ks decay length seen afterwards and Kl inelastic scattering with two charged tracks 
seen in the final state. 

The combination of the described processes was used to fit the experimental dis- 
tribution and seems to describe it well enough. The vertex spatial resolution obtained 
from the Ks decay length distribution was taken into account. The radius at which 
the reconstruction efficiency begins dropping down was found to be 5.0±0.4 cm with 
a 4.8±0.4 cm visible exponential decay length. 

The number of Kl having interactions in the Be tube was found to be 79 ± 18 with 
1355 events representing Kl decaying in flight. 

To select candidates to Kl -^ tt+tt^ events, an additional cut requiring the invariant 
mass of two tracks from a Kl vertex to be in the range of 470-525 MeV was applied. 
The decay radius distribution is presented in Figure 4d together with the flt function 
where all parameters except the number of events are flxed at the values obtained from 



the distribution in Figure 4b. The number of events under the peak drops down to 31 
± 7 and 78 remain from decays in flight. 
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Figure 3: Study of KsKl coupled decays: Invariant mass vs. missing momentum for 
1-st (a) and 2-nd (b) vertex; c. Invariant mass for Ks (shaded) and K^ after Ks 
selection; d. Missing momentum for Ks (shaded) and Ki after Ks selection. 



The peak events after the invariant mass cut are interpreted as regeneration of K^ 
into Ks with its decay into tt+tt^ and are used for the calculation of the regeneration 
cross section. 

For the K^ -^ 7i~^7i~ selection one can apply stronger requirements for these events 
to satisfy full kinematics within detector resolution, i.e. 80 MeV/c< P^ < 140 MeV/c 
and Ks vertex in the Pmis direction. This selection is shown in Figure 4d by the shaded 
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histogram and demonstrates that the peak at the Be tube survives with 20±5 events 
and 35 Kl decays in flight remain, still 10 times more than the expected number of 
CP violating decays. 
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Figure 4: Study of KsKl coupled decays: a. Decay length for Ks'^ b. Decay radius 
for Kl, c. Difference in the angle of Pmis and vertex-beam line for "tube" events 
(histogram), K^ semileptonic decays (points with errors and Ks two pions decays 
(shaded); d. Decay radius for K^s after Mj„„ cut and after Ks selecting cut (shaded). 



The DC material (120 /i thick mylar entrance window, 100-150 fi diameter 50%Cu- 
50%Ti wires, Ar gas) can also contribute to nuclear interaction events. The DC mylar 
entrance window at 2.1 cm radius adds about 10% to the tail of events peaked at the 
Be pipe. The biggest concentration of the DC fleld wires is between 2.5-3.5 cm, where 



the average thickness of Cu+Ti is estimated as 3 x lO^^g/cm^ and drops down with 
radius. This amount of material is neghgible compared to Be and could produce less 
than one nuclear interacting event, but a 10% correction from the mylar window is 
applied to cross section. 

A further correction to the K^ cross sections comes from the A and E hyperon 
decays which can appear in the kinematic range as K^ — > 2tt. The reactions 

^L + P -* A° + 7r+, A° ^ j9 + TT"; 

i^L + j9 ^ SO + 7r+, SO ^ AO + 7, AO ^ p + tt"; 

Ki + n —^ S^ + 7r+, S^ -^ n + tt"; 

i^L + n ^ S+ + TT", S+ -^ n + 7r+ 

were studied. In our selection only two charged pions were detected by the require- 
ment of two prong vertex and dE/dX in the DC corresponding to minimum ionizing 
particles. 

As it was shown by simulation, 8% of such events with a final nucleon nearly at rest 
contribute to the number of regenerated events after the Mj„^ cut with the visible vertex 
radius distribution similar to that from the Ks decay. These events were removed from 
the regeneration candidates in the Kl kinematic region and added to 79 - 31 = 48 Kl 
interaction candidates. 

The following numbers were obtained: 

^^ nucl+reg ly It io, 

Nnucl = 50 ± 17, 
N.,eg = 29 ± 7. 



Events in the range 1.5-3.5 cm from the distribution in Figure 4d were used to 
obtain the angular distribution for Kg after regeneration. Figure 4c shows the projected 
angular distribution of these Ks- Dots with errors show the expected distribution for 
semileptonic decays of K^., normalized to the expected number of these background 
events in the selected range. 

The obtained angular distribution is wider than in the case of coherent regeneration 
which should look like shown shaded distribution for original Ks decaying at the same 
distance. But with this data sample the coherent contamination is expected to be small 
and cannot be extracted. 



4 Cross sections 

For the calculation of the cross sections, the number of the initial Ki passing the Be 
tube multiplied by the reconstruction efficiency e^ec can be found from the distribution 
of Figure 4b using the expression: 

Nk^ ■ erec = No ■ LkJB,h = 61200 ± 2500, 

where A^o is the visible number of pairs at a zero radius (fit value from the fiat region 
in Figure 4b) and L;^^=350 cm is the decay length for the K^. The value i?c/i=0.78 
is the probability for a K^ to have a pair of charged particles after decay. Using this 
expression one can calculate a probability for the Kl to interact at the Be tube as 



-3 
Preg = Nreg/{NK, " e,ec " B^^) = (0.61 ± 0.16) • 10" 



Pnuci = N^uciKNk^ ■ erec) = (0.82 ± 0.28) ■ 10- , 

-3 



5 



where -87^77=0.686 is a branching ratio of Ks decays into a pair of charged pions. A 
ratio (1.12 ± 0.16) for the reconstruction efficiencies of Kl decaying into two pions 
(after regeneration) and into three particles was obtained by simulation and was taken 
into account. 

Using the above probabilities the cross sections can be calculated from the equation 



a- P-^ 



NA-p-t ' 



where A and A^^ are the atomic number and the Avogadro constant, p and t are 
the density and the thickness of the material. After correction for the mylar window 
interactions the following cross section for Be has been obtained: 

crf4 = 53 ± 17 mb. 

For the nuclear cross section (excluding regeneration) one can obtain the inelastic 
cross section into a two particle final state, arising from the A and S production dis- 
cussed earlier: 



T^ei = 70 ± 26 mb. 



To estimate the total cross section, the relative weight of these reactions was found 



to be 0.21 from the NUCRIN package [|16[. Using the CMD-2 experimental crj'^g; value 



and a ratio amei/crtot = 0.52 |T5|, one can estimate a^^ = 641 ± 238 mb. A systematic 
error of this estimation is about 30% and comes from uncertainties in the event selection 
and in the ratio above. 

5 Discussion 

The selection of candidates for Kl -^ tt+tt" events faced two problems. 

First is a background from the dominant Kl decays. Simulation gives a rejection 
factor of only 0.03 for the given DC resolution, and the ratio of the number of events 
within our cuts (away from the beam pipe) to the total number of K^ decays, 35/1355, 
is consistent with this result. However, the expected rate of Kl -^ 2tt decays is 
2 • 10~^, so that only 3 true Kl -^ 27i events could be expected in the sample. The 
observed events are mainly from the semileptonic decays of Kl — *> vr/iz/, and to improve 
the achieved rejection and select Kl -^ tt+tt" events better angular and momentum 
resolution in the DC is needed. 

The effect of tightening resolution may be illustrated by selecting events within 1.5 
standard deviations of the detector resolution. 

In this case 9.5±3.4 of the original events remain under the peak at the pipe and 
only 5 Kl decays in flight remain (2 of them should be real Kl -^ tt+tt^ events). 
DC resolution better by a factor of 2 expected for the rest of our data, could give a 
signal/background ratio about 1 or better. 

A second problem is the relatively high background from nuclear interactions of Kl 
and the regeneration effect. The obtained cross sections can be compared to the data 
available at higher momenta and with theoretical expectations. There are no data for 
the regeneration cross section for Be and for slow kaons. 

In Figure ^ the experimental regeneration cross section is plotted together with 



the theoretical calculations performed in [15|. The obtained cross section is consistent 



with calculations. The calculated regeneration cross section for Cu is also presented in 
the Figure. The comparison of the calculated regeneration cross sections for these two 
different materials shows that at momenta below 200 MeV/c one cannot scale them by a 
simple A^^^ dependence. To verify these theoretical calculations one needs experimental 
measurements for different materials. The experimental angular distribution of the 
regenerated Ks after background subtraction is presented in Figure |p3 together with 
the fit function and the theoretical prediction []T3| and seems to be a little narrower 
than the expectation. 

This result for the total nuclear cross section in Be for 110 MeV/c kaons can be 
compared to the experimental data at higher momenta |]T^ and theoretical calculations 
||T5| . Such comparison is presented in Figure 6. The cross sections extracted from 
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Figure 5: a. Experimental regeneration cross section and theoretical calculations for 
Be and Cu; b. Angular distribution of the regenerated Ks with the best fit function 
(solid line) and theoretical prediction (dots); 



GHEISHA and FLUKA simulation codes are also shown. It is seen that the FLUKA 



code, as well as the calculations from ||15[ give cross sections in good agreement with 
experimental data. The GHEISHA code gives completely wrong absolute values as 
well as momentum dependence. 



After publishing our preliminary results [T^, the regeneration influence was dis- 
cussed for the e'/e measurement planned in KLOE detector [|l^. A possible way to 
remove regeneration events is to require that the final state particles lie in a plane 
(since for regeneration events, recoil nucleons will carry a momentum). It was shown 
that the total regeneration probability in the KLOE drift chamber after an acoplanarity 
cut (factor of 4 rejection) was 10~^ that should be compared with 2 x 10~^ probability 
for the "normal" CP violating Kl —>■ vrvr decay and ~ 10^^ probability for the expected 
direct CP violation decay. 

The regeneration itself does not give any decay asymmetry expected for the direct 
CP violating Kl decay, but the acoplanarity cut as well as other selection cuts applied 
separately to Tr~^TT~ and 7r°7r° final states with a different resolution (DC for the first 
and calorimeter for the second) can cause a systematic asymmetry due to the broad 
angular distribution of the regenerated events. 
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Figure 6: Comparison of the total experimental Kl nuclear interaction cross section 
in Be with the theoretical calculations and simulation by different codes. 



6 Conclusions 



With the CMD-2 detector the coupled KsK^ decays from the have been observed 
for the first time at an e^e~ machine. 

The attempt to select Kl -^ tt+tt^ events emphasized again problems of the Kl -^ 
TTfiu decay mode background as well as a high level of nuclear interactions of neutral 
kaons including Kl into Ks regeneration. 

The regeneration and total nuclear interaction cross sections for low momenta neu- 
tral kaons have been measured. The measured values of the cross sections indicate that 
regeneration will cause an additional background for the CP-violating decays of Kl at 
the 0-factory experiments and should be carefully studied to avoid systematic errors. 

Comparison of the obtained total cross section with the predictions of the simulation 
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codes GHEISHA, FLUKA and calculations performed in |T5| shows that GHEISHA 
results are wrong while two other codes are in good agreement with experimental data. 

The opening of this kinematic region for studies of the neutral and charged kaon 
interactions has been an additional argument for the construction of 0-factories and 
the detector FINUDA ||T^ has been proposed for these studies at the 0-factory in 
Frascati^]. We anticipate that the results obtained from the data now in hand will be 
important in planning of experiments at 0-factories. 

The analysis of new experimental data is in progress and we expect new results on 
the nuclear interaction cross section. 

7 Acknowledgements 

One of the authors (E.P.Solodov) would like to thank R. Baldini for useful discus- 
sions on the nuclear interaction studies and cooperation in this field and A. Michetti 
for simulation of the hyperon production of neutral kaons and theoretical cross section 
calculations. 

This work is supported in part by the US Departament of Energy, US National 
Science Foundation and the International Science Foundation under grants RPTOOO 
and RPT300. 

References 

[1] V.N.Bayer, ZETFP 17 (1973) 446. 

[2] G.A.Aksenov et al.. Preprint BudkerlNP 85-118, Novosibirsk, 1985. 

[3] V.V. Anashin et al.. Preprint BudkerlNP, 84-114, Novosibirsk 1984. 

[4] P. Eberhard, Contribution to the Factory Workshop at UCLA, April, 1990.; 
G.Ghirardi, Proceedings of the Workshop on Physics and Detectors for DA$NE, 
Frascati, April, 1991, p.261. 

[5] J.L. Rosner, I. Dunietz, J. Hauser, Phys. Rev. D35 (1987) 2166. 

[6] A partial list of other major contributions includes: A.N.Skrinsky et at, "Novosi- 
birsk 0-factory project", F.J. Botella, J. Bernabeu and J. Roldan -Blois CP Vi- 
olation Conference, FTUV/89-35,IFIC/89-ll, May 1989.; G. Barbiellini and C. 
Santoni, CERN-EP/89-8 and CERN-PPE/90-124; J.A. Thompson, University of 
Pittsburgh preprint PITT-90-09; contributions by J.A. Thompson, D. Cline, and 
by C. Buchanan and others to the (p Factory Workshop at UCLA, April, 1990.; 
Y. Fukushima, et al,.KEK Preprint 89-159; summary talks by P. Franzini and 

13 



M. Piccolo and others at the Workshop on Physics and Detectors for DA$NE, 
Frascati, April, 1991. 

[7] A.N.Skrinsky, Proceedings of the Workshop on Physics and Detectors for DA<I>NE, 
Frascati, April, 1991, p. 67. 

[8] G. Vignola, Proceedings of the Workshop on Physics and Detectors for DA$NE, 
Frascati, April, 1991, p. 11. 

[9] A.N.Filippov et al.. Proceedings of the XVth International Conference on High 
Energy Accelerators, Hamburg, Germany, 1991, World Scientific, V. II, p. 1145. 

[10] S. Eidelman, E. Solodov, and J. Thompson, Nuclear Physics B (Proceedings Sup- 
plement) 24A (1991) 174. 

[11] E.V. Anashkin et al., IGFA Instrumentation Bulletin 5 (1988), p. 18. 

[12] R.R. Akhmetshin et al., Phys. Lett.B 364(1995) 199. 

[13] R.R. Akhmetshin et al. Preprint BudkerlNP 95-62, Novosibirsk, 1995. 

[14] R.M.Barnett et al., Phys.Rev.D54( 1996)1. 

[15] R. Baldini, A. Michetti, Preprint LNF-96/008, 1996. 

[16] This value was extracted from GERN GEANT package. 

[17] G. A. Sayer and E. F. Beall, Phys. Rev. 169(1968)1045. 

[18] M.Agnello et al. (FINUDA collaboration), FINUDA, proposal LNF-93/021 (IR) 
11 Maggio 1993. 



14 



